Spin-Fluctuation-Driven Orbital Nematic Order in Ru-Oxides: 
Self-Consistent Vertex Correction Analysis for Two-Orbital Model 
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To reveal the origin of the "nematic electronic fluid phase" in Sr3Ru2 07, we apply the self- 



consistent vertex correction analysis to the 



, dj/zVorbital Hubbard model. It is found that the 



Aslamazov-Larkin type vertex correction causes the strong coupling between spin and orbital fluctu- 
ations, which corresponds to the Kugel-Khomskii spin-orbital coupling in the local picture. Due to 
this mechanism, orbital nematic order with C2 symmetry is induced by the magnetic quantum criti- 
cality in multiorbital systems, while this mechanism is ignored in the random-phase-approximation. 
The present study naturally explains the intimate relation between the magnetic quantum criticality 
and the nematic state in Sr3Ru207 and Fe-based superconductors. 
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Recently, emergence of orbital (or quadrupole) or- 
der or orbital fluctuations in multiorbital systems has 
been attracting great attention. In heavy-fermion sys- 
tems, CeBg exhibits non-magnetic quadrupole order [1], 
and the hidden-order phase in URu2Si2 [2, 3] is ex- 
pected as quadrupole or higher-rank multipole order. 
As for d-electron systems, Fe-based superconductors ex- 
hibit "non-magnetic" orthorhombic structure transition 
at T = Ts as well as the nematic order, which indicates 
the occurrence of the ferro-orbital polarization n xz ^ n yz 
[4, 5]. Large softening of shear modulus Cqq above Tg in- 
dicates the existence of strong orbital fluctuations [6-8]. 

The quantum critical phenomenon in bilayer per- 
ovskite Sr3Ru2C>7 is very unique in that both spin and 
charge degrees of freedom are intimately related [9-12]. 
The band structure of Sr3Rii207 is composed of the ti g - 
orbital (d xz , d yz , d xy orbitals) of Ru ions. Both the T- 
linear resistivity [9] and the NMR result \jT x T oc T^ 1 
[11] are observed below 20 K under the magnetic field 
H c 1=3 7.8 Tesla. These non-Fermi liquid behaviors indi- 
cate the emergence of "antiferromagnetic quantum crit- 
icality (AFM-QC)" at H m H c , although no long-range 
magnetic order is observed till 0.1 K [10, 11]. 

Under the critical field H k, H c , moreover, Sr3Ru2 07 
exhibits a novel "non-magnetic nematic electronic fluid 
phase" below 1 K, which is confirmed by the large 
anisotropy of in-plane resistivity [10]. As a possible ori- 
gin, the Pomeranchuk instability of the Fermi surfaces 
(FSs) in the single-band Hubbard model had been stud- 
ied using the renormalization group method [13, 14] and 
the perturbation theory [15]. Also, the orbital polar- 
ized state has been studied based on multiorbital Hub- 
bard models, using the mean-field-level approximation 
(MFA) [16-18]. However, many of these studies do not 
prove the stability of the nematic order against antifcrro- 
spin/orbital order driven by the nesting of the FS. Espe- 
cially, the key question — why the non-magnetic nematic 
order occurs only near the AFM-QC in Si'3Ru2 07 — has 



been still unclear. 

The nematic order in Fe-based superconductor also at- 
tracts increasing attention, and it would offer us useful 
hints in the study of S^R^O^. Since it cannot be ex- 
plained in the MFA analysis, the spin nematic order due 
to the order-by-disorder mechanism had been proposed 
[6, 19]. However, it could not be applied to Si^R^C^ 
since the incommensurate spin fluctuations are realized. 
On the other hand, we have recently revealed that the 
vertex correction (VC), which describes the many-body 
effect beyond the MFA, induces strong ferro-orbital fluc- 
tuations [20, 21]. Then, it is highly required to analyze 
the significance of the VC in Ru-oxides. 

In this paper, we study the origin of nematic phase 
in Sr3Ru207 based on the two-orbital Hubbard Hamilto- 
nian. We utilize the self-consistent VC (SC-VC) method, 
which was recently applied to Fe-based superconductors 
successfully, and reveal that the "orbital Pomeranchuk 
instability" is generally induced near AFM-QC, owing to 
the spin-orbital coupling given by the VC. The present 
study predicts the realization of orbital nematic order 
near the field-induced AFM-QC in Sr 3 Ru 2 7 . 

The bandstructure of Sr3Ru207 is rather complex be- 
cause of the tilting of RuOe octahedron. Therefore, we 
study the simplified (d xz , dj, z )-orbital model to grasp the 
essential mechanism of the orbital nematicity: 



H = 



E 



£fe C fe, J u.,<7 C *W + He 



(1) 



where [i,v = 1,2 represents the d-orbital; 1 = xz and 2 = 
yz. This model describes the FSa and FS/3 of Ru-oxides, 
and it was analyzed in the study of anomalous/spin Hall 
effect [22]. Hereafter, we promise that x, y-axes are along 
the nearest Ru-Ru bond directions. Then, the intra- 
and inter-orbital hoppings are given as ^ = — 2tcosk x , 
= — 2t cos k y , and = At' sin k x sin k y . In this pa- 
per, we put (t,f) = (1,0.1). The FSs for the electron 
filling n = 2 are shown in Fig. 1 (a). He represents 



2 



the multiorbital Coulomb interaction composed of intra 
(inter) orbital interaction U (U') and the exchange inter- 
action J. [21]. Hereafter, we put U = U' + 2 J. 



The divergence of x x 



(0) immediately leads to the 
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FIG. 1: (color online) (a) FSs of the two-orbital model for n — 
2. The colors correspond to 1 = xz (red) and 2 = yz (blue). 
The arrow represents the xz intra-orbital nesting vector. The 
dashed lines represents the deformed FSs for {O x 2_ y 2} < 0. 
(b) Correspondence between (O x 2_ y 2) and d-orbital. (c) otQ, 
o.q and a c for n = 3.3 in the RPA. (d) X s (q) and (e) x"(q) 
obtained by the RPA. 

In the present model, the susceptibility for the charge 
(spin) channel is given by the following 4x4 matrix form 
in the orbital basis: 

r (s) {q) = £ta,cM (g)(1 _ fc( S )£irr,c( S ) (g)) -l ; (2) 

where q = {q,0Ji = 2ttIT), and T c ^ represents the 
Coulomb interaction for the charge (spin) channel com- 
posed of U, V and J given in Refs. [21]. The irreducible 
susceptibility is 



X i ">< s Hq) = X°(q)+X c ^(q), 



(3) 



where Xu', mm '{l) = ~ T E P G im {p + q)G m , v {p) is the 
bare bubble, and the second term is the VC that is ne- 
glected in the RPA. 

In the present discussion, it is convenient to 
consider the quadrupole susceptibilities: Xjil) = 

Ew, m m'O l /xh, mm '(l)Of' m , where &>™ = (l\6 7 \m) 
is the matrix clement of the 7-quadrupole operator. Non- 
zero matrix elements of the quadrupole operators are 
Olt y 2 = ~0 2 x t y 2 = 1 and = 0% = 1 [21], while 
O xz ,yz = in the present Hilbert space. 



ferro-quadrupole order (O x 2_ y 2) 



'■yz 



0, result- 



ing in the "nematic" deformation of the FSs shown in the 
dashed lines in Fig. 1 (a). The correspondence between 
the quadrupole order and d-orbital is depicted in Fig. 1 
(b). The director of the nematicity is along the Ru-Ru 
bond direction, which is consistent with experiments. In 
case of (O xy ) ^ 0, the director is rotated by 45°. 

Here, we introduce the charge (spin) Stoner factor 
otq S \ which is the largest eigenvalue of f c ( s );£ 1IT < c ( s )(q) 
at u\ = 0: Then, the quadrupole (spin) susceptibil- 
ity is enhanced in proportion to (1 — amax) -1 , where 
a.maL = maxq{ag S '}. Since J > in real systems, spin 
fluctuations are always dominant (a„ lax > a^) in the 
RPA as shown in Fig. 1 (c): This figure shows the JJU- 
dependences of q, c q and q,q in the RPA, where U is de- 
termined by the condition ct s Q (U) = 0.97. In the SC-VC 
method, however, the opposite relation < o;^^ can 
be realized even for J/U < 0.1 because of large X c (q). 

First, we perform the RPA calculation for n = 3.3 
and T — 0.05, using 64 x 64 fc-meshcs. The unit of en- 
ergy is t = 1. Figure 1 (d) shows the spin susceptibility 
X s (q) = T,i, m Xu, mm (Q) and the ( e ) quadrupole suscep- 
tibility X c x 2_y2{q) for J/U = 0.06. The Stoner factors are 
"max = 0-97, a c Q = 0.64, and a c = 0.42; see Fig. 1 (c). 
Both x s (q) and x x 2_ y 2(q) have peaks at (q*,q*), where 
q* = 0.317T. Thus, the RPA cannot explain the nematic 
transition that requires the divergence of x^ 2 _ l) 2(0). 

In the next stage, we study the role of VC due to 
the Maki-Thompson (MT) and Aslamazov-Larkin (AL) 
terms shown in Fig. 2 (a) of Ref. [20]. They are 
given by the Ward identity T 1 = S^ FhEX /SG using the 
FLEX self-energy. Moreover, they correspond to the 
first-order mode-coupling corrections to the RPA suscep- 
tibility: The intra- (inter-) bubble correction gives the 
MT (AL) term [23]. In single-orbital models, these VCs 
had been studied by the self-consistent-renormalization 
(SCR) theory [23]. However, significant role of the AL- 
type VC in multiorbital systems had been overlooked un- 
til recently [20]. 

The charge AL term X Ah ' c (q) = A >AL - t ' t (g) + 
X Ah '^(q) is given in eq. (5) of Ref. [20].^ Also, the 
II', mm' component of the spin AL term X AL ' s (q) = 
X AL '^(q) - X AL ^^(q) is given as 

I E E A«'.aM/fo k )HV a % cd (k + Q)V e c figh (-k) 

(q;k) 

+^ b<cd (k + q)V e s ftgh (-k)A'^ m , iCdigh (q; k)] (4) 

where V s < c (q) = t s ' c + f s > c x s ' c (q)T s ' c . A(q; k) is 
the three-point vertex [20], and A' mm , tCdtgh (q;k) = 
A c h,mg,dm'(q;k) + A gdt7nCthm > (q; -k - q) and 

^■mm',ed,ghW ^) = ^ch,mg,dm' (q', k) ~ A g d,mcMm'{q] ~k — 

q). Since ! AL < C (0) ~ £ fc A(3V"'(ife) 2 + V c (k) 2 )A', it 



3 



becomes significant when either x c or x s is large, while 
X Ah ' s is less important unless both \ c and X s are 
large since A" <C A'. For this reason, in Ref. [20], we 
have calculated only X c (q) by putting X s (q) = for 
simplicity. In the present study, we also calculate both 
X c (q) and X s (q) self-consistently. Hereafter, we call 
the former (latter) the SC-VC[c] (SC-VC[all]) method. 
We will show that both methods give similar results. 
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FIG. 2: (color online) (a) xl 2 -y 2 (q) an d (b) X n ^(q) given 
by the SC-VC[c] method for U = 1.75 and J/U = 0.06. (c) 
the charge AL term, (e) Qq and Qo in the SC-VC[c] method 
as function of U under the condition J/U = 0.06. (d) Oq 
and ag as function of J/U in the SC-VC[all] and SC-VC[c] 
methods, under the condition a„ lax = 0.97. 



First, we present the numerical results given by the SC- 
VC[c] method: Figure 2 (a) shows X c x i_ y i{q) obtained 
for n = 3.3, J/U = 0.06 and U = 1.75, in which the 
Stoner factors are a^ ax = 0.93 and = 0.97. In com- 
parison with the RPA, Xx^-y^il) ^ s strongly enhanced 
by the charge AL term, while the results are almost 
unchanged even if MT term is dropped. Figure 2 (b) 
presents the momentum dependence of the charge VC, 
X n L {1(q), which is enhanced on lines q x — 0,2q* be- 
cause of the two-magnon process in Fig. 2 (c), reflect- 
ing the good nesting of the quasi one-dimensional xz- 
orbital bands. (Fig. 2 (c) shows a virtual decay process 
of an orbiton into two magnons with opposite momenta 



±Q.) In the same way, A| 2 22 (q) is enhanced on lines 
q y = 0, 2q* . Thus, both X{ x 11 (q) and X% 2 22(9) gi ye the 
enhancement of x x 2 - y 2 (l) at q = 0. In contrast, Xxy(Q) 
is enhanced only slightly since the three point vertex A 
for O^-quadrupole is much smaller. Figure 2 (d) gives 
the [/-dependence of and ol s q under the constraint 
J/U = 0.06 in the SC-VC[c] method. Since the slope of 
«q is larger than that of a S Q, the critical value (J/U) c 
increases in the strong correlation region. 

Now, we discuss the mechanism of the enhancement 
of x^2_ y 2(0) in more detail: Considering the small inter- 
orbital mixing due to t'(<£. t), we drop Xw'mm'i®) m 
eq. (2) except for I = V = m = m! . Then, eq. (2) is 
simplified to 2 x 2 matrix equation with 

X!i,n(0) Xu, 22 (0) \ 

Xii, 22 (o) x5i,u(o) J ' 
xiT;u(o) \ 
xiT;n(o) J ' 

U 2U'-J\ 

2U' -j u r 




(•5) 



For U' = U — 2J, the charge density susceptibility 
Xn(0) = E;, m Xc' mm (0) and X%^ y M are obtained as 

X c n (0) = 2 x ^ 1 (0)(l + (3C/-5J) X i r 1 i ;ti(0))- 1 (6) 

X^-^(O) = 2xiT;n(0)(l-(C/-5J)xiT.n(0))- 1 (7) 

Therefore, X x 2 -y 2 (®) ^ s enhanced by Coulomb interac- 
tion when J/U < 0.2, while XnW ^ s always suppressed. 
If we drop the spin VC, the spin susceptibility is X s (Q) — 
2xu u(Q)(l - Ux°n ii(Q)) -1 - Therefore, the relation 
a c > a s Q holds for 4i,ii(0) > (1 - 5J/tO _1 x2i,ii(Q) - 
X° n n(0) for J/U < 0.2. Since A AL ' C (0) grows in propor- 
tion to Tx s (Q) [log{x s (Q)} 2 ] at high [low] temperatures 
[20], x x 2 -y 2 (0) is strongly enhanced near the AFM-QC 
under the condition J/U < 0.2. This condition is ex- 
pected to be satisfied in Ru-oxides. 

Figure 2 (e) shows the J/C/-dependence of Uq and a S Q 
given in the SC-VC[all] and SC-VC[c] methods, by ad- 
justing U to satisfy a max = max{af nax , a^ax} = 0.97. In 
both methods, the obtained results are similar since X s 
is less important. The critical value of J/U, at which 
a c Q = a s Q = 0.97 is satisfied, is (J/U) c = 0.095 (0.086) in 
the SC-VC[all] (SC-VC[c]) method. Therefore, the rela- 
tion «□ > «q is satisfied for J/U < (J/U) c ~ 0.1 in the 
present SC-VC method, in highly contrast to the PRA 
result {J/U)f PA = 0; see Fig. 1 (c). 

Figure 3 shows (a) x s (<7) and (b) x%*- y *{q) S iven by 
the SC-VC[all] method for U = 1.35 and J/U = 0.06. 
The Stoner factors are Qq = 0.94 and a|j = 0.97. In 
this method, the shape of X s (q) is slightly changed from 
the RPA result by the momentum dependence of X s (q). 
However, the overall results of spin and quadrupole sus- 
ceptibilities are unchanged by the spin VC. The obtained 
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(a) %(q) 



(b) i c x2 Jq) 




FIG. 3: (color online) (a) X s (q) an d (b) Xx 2 -y 2 (Q) obtained 
by the SC-VC[all] method for U = 1.35 and J/U = 0.06. 



divergent behavior of Xx 2 - y 2 (®) ^ s consistent with the 
"non- magnetic nematic phase" in S^R^Oy near the 
field-induced AFM-QC. Recent ultrasonic measurement 
of the elastic constants peports a large softening of shear 
modulus C s = (Cn - Ci 2 )/2 under H z ~ 7.8 Tesla [24]. 
Since Co 



Cs o « Xl 



j(0) (Cs,o being the lattice 
shear modulus), the obtained development of Xx 2 -y 2 W 
in Fig. 3 (b) is consistent with experiment. 

The role of the AL term can be interpreted as the fol- 
lowing effective nematic interaction within the mean- field 
picture: H ncm - g J2k,k> O x 2 - y 2 (k)O x 2_ y i (fc'), where 
6 x 2_ y 2(k) = Y,i, m , a 0^_ y2 cl tl a c k .rn,cr and g is approxi- 
mately given as (-U + 5J)(Xf hll (q) + X^ 22 (q))/2 X °(q) 
at q = 0. according to eq. (7). In previous studies 
of Sr 3 Ru2 07 based on the single-band model, on the 
other hand, the phenomenological Pomeranchuk inter- 
action H' nsm ~ g'J2k fc' ^fc^fc' had been frequently intro- 
duced, where dk = (cosk x — cos k y ) J2 a c it a Ck , cr - They 
are similar in symmetry, since both dk and O x 2_ y 2 (fc) be- 
long to the same B lg representation of D±h point group. 
Thus, "orbital Pomeranchuk interaction" is driven by the 
AL-VC in multiorbital models. 

It is useful to comment on the difference between the 
present study and the study of the five-orbital model in 
Ref. [20]. In both models, we obtain the development of 
X^2_j,2(0) using the SC-VC method. In the latter model, 
in addition, X xz ( yz ){Q) also develop because of the VC 
and good inter- orbital (xz/yz-xy) nesting. Both ferro- 
and antiferro-quadrupole fluctuations induce the s-wave 
superconductivity without sign reversal (s++-wave state) 
[20]. On the other hand, X xz ( yz ){q) wm be small in the 
*2g-orbital model for Ru-oxides because of ill inter-orbital 
nesting. Thus, the results of this paper would be quali- 
tatively unchanged in the three-orbital model. 

Finally, we discuss the physical meaning of the VC 
in the strong-coupling regime (U 3> t): Due to the 
Kugel-Khomskii type spin-orbital exchange coupling ~ 
(t 2 /U)(s l ■ Sj)(O l x 2_ y 2 ■ 3 x2 __ y2 ), the antiferro-spin order 
induces the ferro-orbital order, and vise versa. Such spin- 
orbital coupling should exist also in the metallic state 
(U ~ t), and it is actually described by the AL-type 
VC. For this reason, cooperative development of spin 
and orbital fluctuations is obtained in the SC-VC anal- 



ysis. [The RPA is insufficient in multiorbital Hubbard 
models in that the spin-orbital coupling is completely ig- 
nored.] Especially, non-magnetic orbital nematic order 
can be realized since the scalar order parameter O x 2_ y 2 
is more stable than the vector order parameter s against 
the quantum and thermal fluctuations. 

In summary, we have studied the origin of non- 
magnetic nematic order in Sr3Ru20y, by applying the 
SC-VC method to the two-orbital Hubbard model. We 
have found that the present model exhibits the orbital 
Pomeranchuk instability near the magnetic quantum crit- 
icality, owing to the spin-orbital coupling described by 
the VC. For J/U < {J/U) c ~ 0.1, the ferro-orbital or- 
der (q;q rs 1) occurs prior to the magnetic transition 
(a S Q ~ 1). The present mechanism gives a natural ex- 
planation for the nematic order in Sr3Ru207 as well as 
Fe-based superconductors near AFM-QC, and it will be 
also realized in various multiorbital models. As the origin 
of the field-induced AFM-QC. both the van-Hove singu- 
larity [26, 27] and the field-suppression of quantum fluc- 
tuation [28] mechanisms had been discussed. 

We note that the renormalization group (RG) method 
is very powerful for the study of VCs in low dimensional 
systems. Recently, we have performed the RG analy- 
sis for the (c? a;z ,(ij /z )-orbital Hubbard model [29], and re- 
vealed that x x 2_ y 2{Q) is critically enhanced by both the 
magnetic and superconducting QCs. 

This study has been supported by Grants-in-Aid for 
Scientific Research from MEXT of Japan. Part of numer- 
ical calculations were performed on the Yukawa Institute 
Computer Facility. 
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